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ABSTRACT 

We study the chemistry in the slow wind during the transition from AGB star to 
pre-planetary nebula (PPN) and planetary nebula (PN). We show that there is a 
very rich chemistry of degradation products created by photoprocessing, driven by 
the gradually hardening radiation field of the central star. Most of these products are, 
however, photodissociated during the PPN phase if the wind is smooth. By contrast, if 
the wind is clumpy, possibly because of dumpiness in the AGB atmosphere, then many 
of these degradation products survive into the PN phase. Thus, chemistry may be used 
to infer the existence of dumpiness in the AGB phase. We identify potential molecular 
tracers, and we note that, in the case of dumpiness, large molecules may survive the 
transport from the stellar atmosphere to the interstellar medium. We compare between 
our model results with observations of three objects at differing evolutionary stages: 
CRL618, NGC 7027 and the Helix nebula (NGC 7293). 

Key words: ISM: globules - ISM: molecules - planetary nebulae: general - planetary 
nebulae: individual (Helix nebula, NGC 7027, CRL618) 



1 INTRODUCTION 

Nearby planetary nebulae are clearly observed to be 
clumpy iO'Dell et al.l 12002') with pe rhaps the best known 
example beina the He lix nebula jMeaburn et alj Il998l : 
lO'Dell fc HandronllTQQ^ . The latter authors estimate that 
around half of the mass of the Helix nebula is contained 
in clumps. The origin of the clumps remains very unclear 
and the re have been s everal models proposed for their gen- 
eration. iDvson et al.l il989l) suggested that the SiO maser 
spots observed in the atmospheres of AGB stars could be 
the first manifestation of clumps later observed in the plan- 
etary nebula. Other models involve clum ps forming much 
later during the planetary nebula phases. Iwilliama (Il999l) 
has shown how density fluctuations at an ionization front 
boundary can produce knots and tails when the ioniz ation 
front switches locally from R-type to D-type (see e.g. IShul 
I1991I for a detailed discussion of ionization front types). In- 
stabilities in the shell of ion ized gas, with the instability 
being cither Raylcigh- Taylor l'Caprio ttilll973l) or 'thin layer' 
( Garcia- Scgura & Franco 1996 ) , have also been proposed as 
mechanisms to form clumps. 

Very recently. iHuggins fc MauronI i2002l) have observed 
the envelopes of NGC 7027, the proto-typical young plane- 
tary nebula, and IRC-1-10216, the standard C-rich AGB star, 
in order to look for knots generated at the AGB phase. Since 
none were observed in these two objects, any envelope knots 



must be very small or else they would be seen in scattered 
light. However, in the innermost regions of NGC 7027, the 
ionized and neutral material surroundi ng the ionized cav- 
ity are observed to be very clumpy and lHuggins fc MauronI 
l|2002) suggest that the knots are formed as the complex 
(and poorly understood) transition from AGB wind to plan- 
etary nebula occurs. The models and explanations for the 
origin of the clumps can thus be separated into two types: 
those in w hich the clumps form early, before the PN is fully 
developed JPvson et alll989l:'Huggins fc Mauron^2002l) and 
those in which the clumps form later during the P N stage 
('Wi Uiamj Il999l : ICarcfa-Segura fc Francol Il996l : ICapriottil 



The chemical history of the clumps will depend on 
where and when they formed since the chemistry is sen- 
sitive to the radiation field, temperature and extinction of 
the clump gas. Molecular observations may thus offer the 
best way to establish the evolutionary history of the clumps 
and thus to differentiate between the two types of model. 
There has been rapid progress in recent years in charac- 
terising the chemistry of the post-AGB phases with many 
molecules being identified in youn g objects such as NGC 
7027 (e.g. iHaseg awa fc Kwokll20oJ and references therein) 
and several species obse rved in evolved PN e such as the 
Helix where, for example. lYoun^^^ (^9 7[l obse rved neu- 
tral carbon in the globules. iBachiller et al.l (1199711 observed 
molecules in clumps in a range of planetary nebulae at dif- 
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ferent evolutionary stages, allowing trends in the molecular 
abundances with age to be identified. 

In this paper, we trace the chemical evolution of model 
clumps, assuming them to have formed early on in the his- 
tory of the PN and compare the results from those expected 
if the clumps formed in situ. The aim is to identify observ- 
able species that can be used as discriminants between the 
two scenarios. 

Our model is described in Section 2, and our results are 
presented, discussed and compared with observations in Sec- 
tion 3. We conclude in Section 4 that molecular tracers do 
exist in PPNe and PNe, and that some large molecules orig- 
inating in the stellar atmosphere may survive the transition 
into the interstellar medium. 



2 PHYSICAL AND CHEMICAL MODEL 

The difference between a clumpy and non-clumpy model is 
that clump interiors will have a higher extinction than their 
surroundings. In the harsh environment of a planetary neb- 
ula this can help to shield molecules so that they survive for 
longer than if they were in the interclump gas. The chem- 
ical difference between clumps that are formed out of the 
planetary nebula gas and those that formed earlier is that 
in the latter case, complex molecules may be shielded and 
preserved long enough to be detectable in clumps. In the 
former case, all such molecules will have been destroyed as 
the extinction o f the g as dropped as the PN evolved. 

iHowe et al.l lll992l) investigated the gas-phase chemistry 
in a carbon-rich AGB wind during its transitio n to a plane- 
tary nebula. This work was followed by a paper jHowe et alJ 
Il994^ investigating the formation of molecules in a dense, 
neu tral globule such as those observed in the Helix neh- 
ula. lHowe et"ai] Jl994 ran equilibrium chemistry models for 
clumps with extinctions between and 2.0, c omparable to 
those of the Helix knots iMeaburn et al.lll998l) . They found 
that molecular abundances in the globule are enhanced with 
C2H and ON at abundances that are possibly detectable in 
carbon-rich globules. ^^^^^^^ 

Here, we extend the work of iHowe et alJ 1^3) in that 
the time-dependent chemistry is calculated, beginning with 
the rich chemistry of a carbon-rich AGB atmosphere. The 
molecular species are assumed to be locked into a clump 
with a density contrast with the surrounding medium. The 
chemistry of the clump is then traced as both it and the sur- 
rounding medium expand away from the star. At the final 
stage of the calculation, the clump and medium have similar 
properties to those of the Helix. As the knots in the Helix 
have a low extinction, the chemistry by this stage is tend- 
ing towards that expected of a diffuse cloud with this initial 
composition. However, at intermediate stages, the rich initial 
chemistry should give rise to marked differences compared 
with that expected from a standard ISM mixture. Some of 
these molecules should be present in detectable quantities 
and if found, would argue for an early clump formation 
model. We follow the physical and chemical evolution of a 
parcel of gas (either clump or interclump material) as it 
moves out from the AGB atmosphere to the protoplanetary 
(PPN) and planetary nebula (PN) phases, ultimately merg- 
ing with the interstellar medium. 

In the AGB atmosphere, the number density and tem- 



perature are on the order of 10^^ cm~^ and 10^ K, respec- 
tively. We assume that the clump /in terclump density ratio 
is about 10 llMauron fc Hugginj200ol find the multiple shells 
in the envelope of IRC-l-10216 to have density contrasts of 
< 10), and that the size of the clump is comparable to a 
stellar radius (about 10^"^ cm). Thus, if the gas:dust ratio 
is similar to that in the interstellar medium, the visual ex- 
tinction associated with a clump in the atmosphere is very 
large. The chemistry in the atmosphere is determined un- 
der thermodynamical equilibrium, and the relative abun- 
dances are fixed as the envelope expands and the density 
falls. At lower densities, however, thermodynamic equilib- 
rium no longer applies and chemistry is driven by two-body 
reactions of ion-molecule and other types, as in conventional 
interstellar chemistry. 

Since no chemical evolution occurs in the initial expan- 
sion, we begin our computation of the evolution at a later 
stage (time t — ti = 100 yr) at which the densities are low 
enough to be dominated by two-body reactions initiated by 
cosmic ray ionisation and by photoprocesses driven by the 
radiation field of the central star and by the ambient inter- 
stellar medium. 

For simplicity, we use a power-law description of the 
physical development. We do not attempt to model the ini- 
tiation and evolution of a clump. In a steady outflow at con- 
stant velocity, the density should fall off as t~^. However, the 
density fall off in a clump may be slower. A t~^^^ dependence 
gives a reasonable fit to typical densi ties through the PPN 
and PN stages (as observ ed by, e.g. iMartin-Pintado et alJ 
ll995l:lMeaburn et al.llT99ll . Therefore we adopt 
-3/2 



n{t) = 10 



(1) 



for times t > ti, where ti = 100 yr. On the assumption of 
constant clump mass throughout the expansion, this gives 

t \ 1/2 
— cm. 

The visual extinction associated with a clump must therefore 
vary as , and we write 



r{t) = 10^ 



(2) 



10' 



(3) 



The distance of the clump from the star is d, which in a 
steady flow must increase linearly with time; 



.,„.io-(l) 



(4) 



We adopt the following dependence of temperature on time, 
as it gives values for the PN clu mp temperatures comparable 
to the typical measured values jMartin-Pintado et aljir995l : 
iMeabu rn et al. 1998): 



m = 316 ( - 



-1/4 



K. 



(5) 



The radiation field experienced by the parcel of gas is ini- 
tially dominated by the contribution from the central star, 
and later by that of the interstellar radiation fiel d. We adopt 
the fo llowing expression, similar to that used bv lHowe et alJ 
il994l) . for the radiation field intensity x, 

(6) 

+ Xism; D < 10"^ cm, (7) 



X = 0; D< 10"° cm, 
-2 



100x..(^)' 
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where we assume that Vism g ives photorates as in the UMIST 
rate file iMillar et al.lll99^ and discussed further below). 
There is some uncertainty in this, since the radiation field 
may be harder than that of the interstellar medium. We 
use the canonical cosmic ray ionisation rate of Xism = 1.3 x 
10-1^ s-i. 

The evolution of the interclump medium can be followed 
in a very similar way to that of a clump. In this single-point 
calculation, the initial composition in the AGB atmosphere 
is the same as for a clump but the parcel of gas is assumed to 
have an initial density lower by a factor of 10 compared with 
the clump. The density falls off with the square of the dis- 
tance from the star, as appropriate for a steady spherically 
symmetric flow. 

n{t) = 10^ ' cm-^ (8) 

A parcel of gas will then experience an extinction that varies 
as 




Since the clump and interclump medim will initally have the 
the same temperature, it is assumed for simplicity that the 
interclump medium temperature is the same as that of the 
clump, and therefore varies as given by Equation 5. This 
is unlikely to be true in the later stages of the PN evolu- 
tion when the interclump medium will be hotter than the 
clumps but this will only accelerate the destruction of those 
molecules that survived the drop in extinction. In the pe- 
riod of interest, the range of temperatures involved does not 
influence the chemistry very strongly. 

The chemical model u sed is similar to the one employed 
m IViti fc WiUiamsl l^^9^. However, we have made consid- 
erable changes to the chemistry. To all the basic species (O, 
C, CO etc.) we have added all the species (mainly hydrocar- 
bons) that are thought to be abundant at the initial stage 
of the PPN; this chemistry also includes species that have 
an unpaired electron located on an internal or terminal car- 
bon atom. For most of these additions, the UMIST rate file 
for gas phase reac tions, involving several hundred species 
l|Millar et al.lll99/t l. and which is normally used routinely in 
our models, did not have a sufficient network of reactions. 
Therefore the formation and destruction routes were taken 
from Frcnklach & Fcisclson (1989). The reaction rates were 
fitted by using different formulae from the UMIST ones. 
To account for that, we have modified our code to in clude 
new rate coefficients calculations taken from lCaul (|200Jl . For 
some of those reactions, no rate coefficients are available and 
[Frenklach & Fcieolson (1989.) determine them from the rate 
coefficients of their inverse reactions: we have employed the 
same method here. We have not included three-body reac- 
tions. We have also excluded reactions which use an 'unspec- 
ified' species as a coUisional body causing de-excitation or 
dissociation. These may be of some r elevance and should be 
included in further studies. Note that lFrenklach fc Feieelsonl 
l|l989() also extrapolated their rate coefficients to tempera- 
tures very different from those at which they were measured. 



Time 




-12 ' ' ' ' 

100 1000 lO" 

Figure 1. Clump (solid lines) and interclump (dashed lines) 
abundances for CO (uppermost solid line) and CN. The time 
is measured in years from the end of the AGB phase. The frac- 
tional abundance is with respect to the total number of hydrogen 
atoms. The object becomes a PPN after a few hundred years and 
by about 1000 years the object is a PN 

3 RESULTS AND DISCUSSION 
3.1 Predictions 

We present in Figures 1-4 results of the calculations of the 
evolution of molecular fractional abundances in clump and 
interclump gas, for several illustrative species. In Table 1 we 
give the computed fractional abundances of potentially ob- 
servable species in the clump, and also the clump:interclump 
ratio of abundances, for several epochs between about 2 and 
10 thousand years, i.e. the period of transition from PPN to 
PN phases. 

The figures show generally the same qualitative be- 
haviour. A comparison of clump and interclump abundances 
shows that molecules are photodissociated earlier in the in- 
terclump than in the clump gas, as expected. Thus molecules 
survive longer in the clumps, and so the ratio of clump to 
interclump abundances varies strongly. For example, the CN 
abundance ratio varies from 0.26 at 2550 yr to 270 at 10050 
yr, a factor of 1000. 

Figure 1 shows the clump and interclump evolution 
of CO and CN fractional abundances. These represent 
molecules that have high and low initial abundances, re- 
spectively, i.e. in the stellar atmosphere of the AGB star. 
The abundance of CO remains high in both clump and in- 
terclump gas until photodissociation begins to play a role as 
the gas density and extinction decline during the expansion, 
with CO in the higher density clump gas surviving until 
later times. On the other hand, CN is a molecule formed as 
a product of photodissociation of larger species, here called a 
degradation product. Since photodissociation is faster in the 
interclump gas, the production of CN is also faster there ini- 
tially, so CN rises faster in the interclump gas than in the 
clump, but also declines earlier. 

Similar effects are observed for other species. The simple 
carbon chains C4H and C5H are both degradation products 
of larger species, and the time evolution illustrated in Fig- 
ure 2 shows that their abundances peak in the interclump 
gas earlier (by several thousand years) than in the clump 
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Species X(clump) x (fnttTump) ^(clump) x(fnttTump) ^(clump) 
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Table 1. Comparison of fractional abundances of potentially observable species at different times. The density, temperature and extinction 
at these times can be readily calculated from the equations in Section 2 



gas. Figure 3 shows the time evolution of C2H2 (which is 
abundant in the stellar atmosphere) and C4H2, a degrada- 
tion product of larger species. Their behaviours are as for 
CO and CN. Both molecules ultimately decline as the gas 
enters the interstellar medium. 

Molecules with rather slower photodissociation rates 
survive longer. In Figure 4, benzene has a timescale for 
survival longer than the ten thousand year period consid- 
ered here. This raises the possibility that some large and 
relatively photo-resistant species may survive the transition 
into the interstellar medium and may establish a popula- 



tion there (Williams 2003). Note however that our benzene 
chemistry may be incomplete; in the chemistry adopted here 
its formation and destruction are mainly dominated by the 
neutral-neutral reaction of two C3H3 molecules and its re- 
verse. 

Table 1 illustrates the effects of the different photodisso- 
ciation timescales in the expanding gas, and gives snapshot 
values of the clump :interclump abundance ratio at 2550, 
6300, and 10050 yr, for species that are potentially observ- 
able. The earliest epoch is approximately when the inter- 
clump abundances are at a maximum, while the values in 
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Figure 2. Clump (solid lines) and interclump (dashed lines) 
abundances for C5H (uppermost solid line) and C4H. The time 
is measured in years and the fractional abundance is with respect 
to the total number of hydrogen atoms. 
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Figure 3. Clump (solid lines) and interclump (dashed lines) 
abundances for C2H2 (uppermost solid line) and C4H2. The time 
is measured in years and the fractional abundance is with respect 
to the total number of hydrogen atoms. 

the clump of species formed by degradation are low. For such 
species, then the ratio is low. 

At epoch 6300 yr, many of the degradation products 
have declined in abundance in the interclump gas, while 
their abundances in the clump gas are growing, so the ratio is 
generally large. At the latest epoch shown in Table 1, then 
most of the degradation products have declined in abun- 
dance in the clump gas, and have declined even further in 
the interclump gas, so the ratio is in many cases again large. 
However, by this stage the abundances in the clump gas 
may have declined so far as to make many of the species 
undetectable. If we take a fractional abundance of 10~* to 
represent a characteristic detection limit at 10000 yr, where 
the clump extinction is ~ 1 magnitude, then about 70% of 
the species listed are possibly not detectable at this epoch, 
whereeis only about 30 percent are not detectable at 6500 yr. 
Thus, there is a very significant change in the detectable 
chemistry of the clump gas the PN ages. 



Time 
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Figure 4. Clump (solid lines) and interclump (dashed lines) 
abundances for CeH (uppermost solid line) and benzene. The 
time is measured in years and the fractional abundance is with 
respect to the total number of hydrogen atoms. 







Clump 


Interclump 


X > 10-^ 


at 2550 yr 


CO, CS, HON, C3H, 
HC3N, C4H5, SiS 


CO, CN, SiC, C2, 

C2H, SiC2, C4H3 


X > 10-'' 


at 6300 yr 


CO, CS, HON, C3H, 
HC3N, C4H5, SiS 


C4H4 


X > 10-* 


at 10050 yr 


CS, C2, C4H4, C4H5 


none 



Table 2. List of important potential tracers of the clump and 
interclump gas at the three epochs 

The column density of any species is made up of con- 
tributions from both clump and interclump gas. At the ear- 
liest PPN epoch illustrated. Table 1 shows that chemistry 
is largely determined by the interclump gas, while at later 
epochs it is determined in the clumps. Table 2 lists molecules 
that are likely to be the most important tracers of either 
clump or interclump gas, at the three epochs chosen. 

3.2 Comparison with observations 

Here the results will be compared very generally with ob- 
servations of three objects: CRL618, NGC7027 and the He- 
lix nebula. These objects span progressively advanced evo- 
lutionary stages with CRL618 being a PPN, NGC 7027 a 
young PN and the Helix an evolved PN. As discussed im- 
mediately below, comparison with CRL618 is made difficult 
by the complex geometry of the source. The comparisons 
with NGC 7027 and the Helix are more straight-forward. 

Carbon-rich protoplanetary nebulae such as CRL618 
and AFGL2688 have been the subject of several molecular 
studies in recent years. In particular, CRL618 has been used 
as a typical example of an asymptotic giant branch (AGB) 
star in the transition evolving toward the planetary nebula 
stage CHerpin fc C crnicharo 2000,; Cernicharo et al. 2001a,b; 
Woods et al. 2002). Here we attempt a qualitative compar- 
ison of our models with some observations and theoretical 
work on this object. iHerpin &: Cernichard i2000|) show that 
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to reproduce the observed emission of CO, HCN, HNC, H2O, 
OH and O, a 3-component geometry is needed (see their fig- 
ure 3), with a central torus at high density (~ 5 x 10^ cm"'') 
and a gradient of temperatures ranging from 250 to 1000 
K; an extended AGB remnant envelope with a density of 

5 X 10^ cm~^ and temperature of lOOK; and finally the lobes, 
where high velocity gas is emitted from, with a density of 
10^ cm~^ and a two-component temperature of 200 and 1000 
K. In follow up papers, Cernicharo et al. (2001a, b) observed 
methylpolyynes, small hydrocarbons and benzene along the 
same line of sight and conclude that the emission of these 
species must come from the central torus region. 

The scenario we model is different from the 3- 
component geometry but we can attempt a comparison 
with the species emitted from the AGB remnant compo- 
nent where the temperature is close to the one derived 
in our model (see Table 1), although the density is much 
higher. We find that our best match with the PPN CRL618 
occurs at early times ~ few hundred years. At this time 
the density is relatively high: our CO column density is 

6 X 10^* cm~^, while Herpin & Cernicharo (2000) derive 
a value of 7 x 10^* cm"^ and our HCN is 2.5 x 10^'' cm'^ 
only slightly over a factor of 4 less than the observed one; 
our HNC is, however, over 3 orders of magnitude underabun- 
dant with respect to observations. Given the differences in 
the physical parameters adopted, any further comparison 
would be of limited value. 

It should also be noted that following the ISO discovery 
of benzene in CRL618 (Cernicharo et al. 2001b), Woods et 
al. (2002) presented a chemical model of this proto-planetary 
nebula. They show that the physical conditions of CRL618 
are such to encourage an efficient production of benzene. 
However, similarly to Cernicharo et al. (2001b) results, they 
model the dense inner torus and therefore we are unable to 
make a direct comparison with their chemical model. 

The young PN NGC 7027 is molec ule rich with a large 
number of species identified within it. iHaseeawa fc Kwokl 
l|200l|) have observed several species and collated the results 
of previous work. In their table 3, they list measured frac- 
tional abundances. If the emission that iHasegawa fc Kwokl 
i200ll) detected is regarded as originating from a large col- 
lection of clumps then their results can be directly compared 
with ours in Table 1. Adopting an age of around 1000 yr for 
the nebula and comparing with our results gives the fol- 
lowing: CN is predicted to have a fractional abundance of 
9x 10"^" and Haseeawa fc Kwok t2001.) measure 1.7 x 10"*; 
HCN is 6 X 10"^ (predicted) verses > 1.2 x 10"^ (observed); 
and C2H is 6 X 10"^ vs 1.1 x 10"*. Given the simplicity 
of the model and the inexactness of the comparison to be 
within around one order of magnitude of the observational 
results for several species is very encouraging (the exact age 
of the nebula or the filling factor of the clumps can easily 
yield uncertainties of an order of magnitude). 

To compare the results for the later stages with those 
of the He lix, we can use the last columns of Table 1 with the 
results of iBachiller et all lll997f) . The model results for this 
evolved PN do not agree nearly so well with the observations 
as for NGC7027 and CRL618. Some results are not incon- 
sistent with the observations: the model predictions do not 
exceed the upper limits established for SiO and CS, being an 
order of magnitude less; the predicted HC3N abundance is 
several orders of magnitude less than the upper limit set by 



IBachiller et alJ lll997l) . However the models severely under- 
predict the abundances of HCN, HNC, CN and HCO+. This 
indicates that these molecules whose earlier clump abun- 
dances were comparable to the observed levels survive for 
longer than our simple model predicts. 



4 CONCLUSIONS 

The chemistry of the gas in the slow wind ejected from a 
star in the AGB phase is strongly affected by photoprocesses 
driven by the hardening radiation field of the central star, 
as the system is transformed into a PPN and then to a PN. 
If the slow wind density distribution is smoothly monotonic, 
then the degradation product molecules in the gas are con- 
fined to the PPN stage. However, if the slow wind is clumpy, 
perhaps because of dumpiness in the precursor AGB at- 
mosphere, then molecules in clumps are protected against 
photodissociation until significantly later times, and may be 
present in the PN phase of evolution. We have presented lists 
of molecules that may represent tracers of such dumpiness. 
Benzene survives the transition into the PN phase, suggest- 
ing that some long-lived molecules may last long enough to 
be injected into the general interstellar medium. A compar- 
ison with observations suggests that this model can, at the 
PPN and young PN stage, correctly predict the abundance 
of several molecular species to within an order of magni- 
tude. At later stages, the model severely underestimates the 
abundance of some short-chain molecules. Thus the shield- 
ing of molecules in clumps may be even more effective than 
as described in this simple model. 

The original aim of this work was to chemically identify 
when clumps form in PNe. The comparison of the results 
in the early stages with the data from CRL618 and NGC 
7027 supports the hypothesis that the clumps form at an 
early stage; if they formed later out of low extinction gas, 
the abundances of the molecules would be significantly de- 
pleted compared with the values observed. As noted above, 
the evidence of the continuing presence of molecules in an 
evolved PN like the Helix indicates that the protective abil- 
ity of the clumps persists until late stages. Future work in 
this field should include a more realistic physical model for 
the different post- AGB phases to account for the survival of 
molecules in the clumps. It should also be possible to model 
individual sources and thus carry out a much more detailed 
comparison with observations both from existing facilities 
and also from forthcoming facilities such as ALMA. 
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